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Abstract. We consider nonequilibrium critical dynamics of the two-dimensional Ising 
model for which the initial state is prepared by switching on random fields with zero mean and 
variance H . In the initial state there is no magnetic order but the clusters of parallel spins 
have a percolation transition for small enough H . Using heath-bath dynamics we measure the 
relaxation of the magnetization which shows a reentrance in time. Due to cluster dissolution in 
the early time regime there is a decrease of the magnetization, followed by an increase due to 
nonequilibrium domain growth which itself turns to a decrease due to equilibrium relaxation. 
The power law decay of the nonequilibrium autocorrelation function is not influenced by the 
percolation in the initial state. 



Introduction. - In nonequilibrium critical dynamics [1-4] the system under consideration 
is prepared in some initial state from which it is quenched to the critical temperature, T c , and 
then let to evolve in time according to the given dynamical rules. Generally one is interested in 
the relaxation of the magnetization, m(t) = (a(t)), and in the behavior of the autocorrelation 
function, G(t,s) — (a(t)a(s)). Here cr(t) denotes the operator of the order parameter and s 
and t are the waiting time and the observation time, respectively. 

In most of the studied cases the initial state is of two kinds: it is either the (completely) 
ordered one, or the (completely) disordered one. Starting from the ordered state with the 
initial temperature Ti — 0, the critical relaxation process involves only equilibrium critical 
exponents. For example the decay of the magnetization is asymptotically given by 

m(t) ~ r*/*, Tj = , (1) 

where x = (3/v is the anomalous dimension of the magnetization and z is the dynamical 
exponent. The decay of the autocorrelation function follows the rule: 

G(t,s)~(t-s)- 2x ^g(t/s), T^O, (2) 

Where the scaling function g(y) ~ y x / z for y ^> 1 and tends to a constant otherwise. This 
means that G(t,s) decreases as t~ x l z for s -C t, similar to the magnetization in Eq. (JJ, 
whereas for t/s = 0(1) the decay is the same as for the equilibrium autocorrelation function. 
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If on the contrary the initial state is the completely disordered one with T, — > oo, new type 
of nonequilibrium exponents appear in the relaxation process. Starting with a small initial 
value, rrii, the magnetization behaves as [5] 

m(t) ~ rriit 9 , t <U, Tj = oo . (3) 

Here is the initial slip exponent, which is zero in mean- field theory but generally positive 
in realistic systems. Consequently in the initial period, t < ti, the order in the system is 
increasing due to critical fluctuations. The limiting time-scale ti can be obtained from the 
condition that for t > U the relaxation goes over to the decay of Eq. QJ, leading to the 
estimate: 

ti ~ m i z ^ x ' , Xi — 9z + x . (4) 

Here Xi is called the anomalous dimension of the initial magnetization. Note that ti tends 
to infinity as m, goes to zero. Another new feature of nonequilibrium dynamics starting 
from a disordered initial state is given by the fact that the autocorrelation function, which is 
measured in the state with m, = 0, is non-stationary: G{t, s) depends both on t and s and 
not only on the time difference t — s. In the limit t^> s the decay is given by [6]: 

G{t,s) ~ t~ x/z , i>s, T = oo, (5) 

and the autocorrelation exponent A satisfies the scaling relation [7] \ = d — 9z = d— Xi+x. 

At this point we can note a formal analogy between nonequilibrium critical dynamics and 
static critical behavior in semi-infinite systems [8-11]. In the former problem translational 
invariance is broken in time due to the quench, whereas in the latter phenomenon translational 
invariance is broken in space in the direction perpendicular to the surface of the system. The 
initial state in dynamics with = (Tj = oo) corresponds to a semi-infinite system with 
fixed (open) boundary conditions. In this respect the exponent Xi is analogous to the surface 
magnetization exponent x\ at an ordinary surface transition. 

In nonequilibrium critical dynamics there are a few examples in which the initial state 
contains nontrivial power-law correlations. Here we mention relaxation in the 2d XY- model 
[12, 13] when the initial state is prepared at a temperature T\ > and the quench is made 
to a temperature T2 > T\ with both temperatures in the critical phase, i.e. T2 < Tkt, where 
Tkt is the Kosterlitz-Thouless temperature. According to spin-wave theory for Ti <C Tkt 
the autocorrelation function in Eq. (J2J is modified by replacing x by x(Ti) — x(T{), where 
x(T) is the temperature dependent anomalous dimension. This theoretical prediction has been 
checked numerically [14]. In another work the nonequilibrium dynamics of the d-dimensional 
spherical model is studied starting from an initial state which has long-range correlations [15]. 
Depending on d and the correlation exponent a, different types of critical ageing behavior are 
found. In both examples described above there is quasi-long-range order in the initial state 
which is characterized by a fractal dimension d\ larger than the fractal dimension dj = d — x 
of the critical state. 

In the present paper we consider a type of initial states in which the disordering effect is 
not due to temperature, but is due to switching on random fields of zero mean and variance H. 
In the limiting case of strong fields, H ^> J, J being the coupling of the system, we recover the 
completely disordered initial state with Tj = 00. In the other limiting case, H — > 0, the initial 
state is completely ordered as for Ti = 0. We are interested in the behavior of nonequilibrium 
dynamics for intermediate values of A = H/J = 0(1). To be concrete we consider the 
two-dimensional Ising model for which random fields of any strength destroy the magnetic 
long-range order [16]. However the structure of the initial state on the level of clusters of 
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parallel spins shows interesting features [17]. For strong enough fields, A > A perc , these 
clusters have a finite extent, whereas for weaker fields, below A perc , there is a percolation 
phenomenon: there are giant clusters for both spin directions which percolate the sample. 
According to numerical results [17, 18] these clusters are isomorph with those of short range 
percolation having the same fractal dimension [19] d P j erc = 91/48. 

In this paper we want to clarify the effect of the change in the initial state topology on the 
properties of the nonequilibrium dynamics. In particular we want to see how the percolation 
of parallel spins is manifested in the magnetization relaxation. 

Numerical results. - The initial state used in our nonequilibrium relaxation measure- 
ments is generated by switching on random fields, so that the Hamiltonian of the system is 
defined by: 

H = -J22 a i a j -'^2hi<Ji , (6) 

(ij) ' 

where <r, = ±1 is an Ising spin located at site i of a square lattice. The nearest- neighbor 
coupling, J > 0, is ferromagnetic whereas the hi random fields are taken from a symmetric 
Gaussian distribution: 

(7) 

which has a variance H. The ground state of the system for a given realization of the random 
fields is exactly calculated by a combinatorial optimization method [20] which works in strongly 
polynomial time. With help of this effective algorithm systems as large as L — 500 are treated 
numerically and averaging is performed over at least 80000 realizations. 

As we have already mentioned in the introduction the magnetic correlation function in 
the initial state, (<7j<Xj), is short-ranged, so that the magnetic correlation length £ m is finite. 
However, considering the topology of clusters of parallel spins other length-scales can be 
defined, which can diverge by varying the parameter A = H/ J. One of these lengths is 
the so-called breaking- up length, lb, which measures the typical fluctuations of the interface 
separating the different spin orientations. In other words in a finite system with L < Ij, the 
ground state is homogeneous, but for L > lb both spin orientations are present in the ground 
state and one can use a coarse-grained description with effective cells of size ~ lb- For weak 
random fields the breaking-up length asymptotically behaves as [21] 

l b ~ cxp(^/A 2 ) , (8) 

thus it is divergent as A — > 0. The other length appearing in the geometry of the initial state 
is the percolation correlation length, £, perci which is the measure of the typical size of the 
largest clusters in the system. £ perc is divergent below the percolation transition point, which 
is estimated as [17] A perc ss 1.65. 

In the numerical calculations we considered the disorder in such a range for which the 
coarse-grained description is applicable, i.e. the breaking-up length satisfies the relation: 
£,L ^> Zfc(A). In this way we made the calculations for 1.25 < A < 4. For example, for 
A = 1.4, deeply in the percolation regime, the breaking-up length is h ~ 25 [17]. 

Having prepared the system in the initial state with a small magnetization m^, we quench 
it to the critical point and study the nonequilibrium relaxation process, thereby using the 
standard heat-bath algorithm. As shown in Figure^the magnetization displays an intriguing 
reentrance in time for not too large values of A: immediately after the quench the magnetiza- 
tion first decreases (regime 1), reaches a minimum at time t min (the value of which depends 
on A and m;), increases in the regime 2 with t min < t < t max , before decreasing again in the 
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Fig. 1 - Typical temporal evolution of the magnetization in a critical relaxation measurement when 
starting from a RFIM ground state. In order to make the three different regimes and the reentrance 
in time clearly visible, we show here data obtained for a large initial magnetization rru = 0.25. The 
variance of the field is A = 1.7 und the system contains 200 x 200 spins. 

large time limit with t > t max (regime 3). This non-monotoneous behaviour of the magneti- 
zation is understood by looking at the snapshots of Figure At t = the system is formed 
by homogeneous cells wherein the spins have the same orientation. The typical linear size of 
these cells is h. Immediately after the quench these compact cells are dissolved which leads 
to a decrease of the magnetization (see the snapshot for t = 5). As this process is exactly 
the same as when starting from a totally ordered initial state (corresponding to one large, the 
system filling, homogeneous cell) the decay of the magnetization follows the power law QJ. 
This dissolution stops when the total magnetic moment of the cells is of the order 0(1), i.e. 
when the magnetization has been reduced by a factor if where d — 2 is the dimensionality of 
the system. This yields the prediction 

(9) 

and therefore ^ 

lni mm ~ liil b ~ —j. (10) 

After the cells have been dissolved, the usual nonequilibrium domain growth sets in yielding 
an increasing magnetization in the regime 2. This increase of order (see the snapshot in Figure 
121 at t — 150) is due to the critical fluctuations and is therefore similar to that observed when 
starting from a completely disordered initial state with a small magnetization. We therefore 
expect that the temporal evolution of the magnetization in the regime 2 is also governed by 
the initial slip exponent 9, see Eq. (J3J, and that the domain growth continues up to a time 
tmax with In t m ax ~ In t m in + In ti where ti is given in Eq. (j^J . For longer times (regime 3) 
the relaxation again goes over to the well-known decay QJ. 

We have studied this intriguing behaviour of the magnetization in a systematic way by 
varying the variance of the random fields over a large range. Some of our main findings 
are summarized in Figure |3| where we present results obtained for the initial magnetization 
rrij = 0.0390625. We checked that our conclusions do not depend on the chosen value of m, 
by simulating systems with other values of the initial magnetization. 
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Fig. 2 - Snapshots of the spin configuration of a system with 128 x 128 spins and = 0.0390625 at 
three different times: ground state at t=0 with A = 1.7, close to the minimum at t — 5, and in the 
increasing part at t = 150. 



Figure shows the temporal evolution of the magnetization at short times for various 
values of the variance A. Starting with the smallest value one observes that by increasing 
A the minimum gets shallower and its position is shifted to smaller t. At a threshold value 
At « 2.0 (the exact value of A t depends on the value of m^) the minimum disappears and 
the magnetization increases monotonically at early times. As mentioned before, the existence 
of a minimum follows from the competition of two processes, cluster dissolution and domain 
growth. When A increases, the breaking-up length (and therefore also the effective length 
of the cells) decreases, making the cluster dissolution less and less effective, up to the point 
where domain growth is the only relevant mechanism, yielding an increasing magnetization 
immediately after the quench. 

Figure gives a closer look at the location of the minimum. Plotting hit min as a func- 




Fig. 3 - (a) Magnetization vs time for various values of the variance A: 1.4, 1.5, 1.6, 1.7, 1.8, 2.0, 
2.2, 2.6, 3.0 (from bottom to top), (b) Logarithm of t m i n as function of 1/A 2 for 1.25 < A < 1.95. 
For variances A > A c m 1.65 the data follow a straight line with slope 9.55. Deviations appear for 
A < A c . Error bars are smaller than the size of the symbols, (c) Magnetization vs time in the regime 
2 for A = 1.4, 1.6, 1.8, 2.0, 2.2, 2.6, 3.0 (full lines from bottom to top) in a log-log plot. The dashed 
line is obtained when one starts from an infinite temperature initial state. All curves yield the same 
slope for large t. 
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Fig. 4 - Temporal evolution of the autocorrelation function G(t, s = 0) for different values of A: 1.4, 
1.6, 1.8, 2.0, 4.0 (from top to bottom). The grey line results from an infinite temperature initial state. 
The inset shows the same quantity for A = 1.4 and systems with different linear extend L: 128, 200, 
500 (from bottom to top). 



tion of 1/A 2 , we see that the points with A £ [1.65, 1.95] lie on a straight line, in complete 
agreement with the theoretically expected behaviour (|10|) . For values of A < 1.65, however, 
deviations from the straight line behaviour are observed. Recalling that the percolation tran- 
sition takes place at A c 1.65, we interpret these deviations as effects on the properties of the 
noncquilibrium relaxation produced by the underlying percolation transition. The percolation 
point thus introduces a finite time and a finite length scale which are the values of t m i n and 
h at A c . 

Finally, Figure shows the time dependence of the magnetization in the intermediary 
regime 2. For any value of A we end up with a power-law increase with an exponent whose 
value is fully compatible with the value 9 s» 0.19 of the initial-slip exponent of the two- 
dimensional Ising model, see (@J and the dashed line in Figure It may, however, take quite 
a long time before reaching this simple power-law regime, as the crossover time increases for 
decreasing A. 

We end the presentation of our numerical data with a discussion of the autocorrelation 
function. As already mentioned, see Eq. JSJ, the autocorrelation follows a simple power law 
in the long time limit when quenched from infinite temperature to the critical point, with an 
exponent X/ z = dj z — 9. As shown in FigureOt relaxation measurements yield in the regime 
2 a common, A independent, value of 9. We therefore expect to regain for the autocorrelation 
the same value for X/ z for any value of A, even so the non-trivial initial state may have a long 
lasting impact on the spin configurations. This is indeed what we observe, as shown in Figure 
0] Plotting lnG(i,0) as a function of hit yields at late times straight lines with a common 
slope X/z w 0.73. Corrections at early times are more important for smaller values of A, 
i.e. for larger effective cell sizes in the initial state, as expected. We finally mention that for 
small values of A strong finite-size effects appear, see the inset of Figure 0J which forced us 
to simulate systems containing as many as 500 x 500 spins. 
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Discussion. - We have reported in this work an intriguing reentrance in time in critical 
relaxation measurements which start from a non-trivial initial state given by a ground state 
of the RFIM model. Competition between two different mechanism, the cluster dissolution 
of the compact cells and the usual domain growth, is responsible for this novel feature in 
nonequilibrium critical dynamics. Note that by varying the strength of the random held, A, 
and the value of the initial magnetization, m,, one can tune the borders of the different regimes 
in Figure n so that both t m i n and U can be macroscopic. Interestingly, the percolation point 
in the initial states, given by a critical variance of the random holds, has an impact on the 
nonequilibrium dynamics by introducing new finite time and length scales. 

Finally, we note that nonequilibrium relaxation based on random fields can be applied 
for other systems, too. For example for the 3d Ising model the initial state is either ordered 
(for weak random fields) or disordered (for strong random fields) [22]. Most interesting is, 
however, the borderline situation when at the critical value of the strength of the random 
fields the initial state is critical and described by a random fixed point. 

* * * 
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